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We identify a novel route to the birth of a strange nonchaotic attractor (SNA) in a quasiperiod- 
ically forced electronic circuit with a nonsinusoidal (square wave) force as one of the quasiperiodic 
forces through numerical and experimental studies. We find that bubbles appear in the strands of 
the quasiperiodic attractor due to the instability induced by the additional square wave type force. 
The bubbles then enlarge and get increasingly wrinkled as a function of the control parameter. 
Finally, the bubbles get extremely wrinkled (while the remaining parts of the strands of the torus 
remain largely unaffected) resulting in the birth of the SNA which we term as the bubbling route 
to SNA. We characterize and confirm this birth from both experimental and numerical data by 
maximal Lyapunov exponents and their variance, Poincare maps, Fourier amplitude spectra and 
spectral distribution function. We also strongly confirm the birth of SNA via the bubbling route by 
the distribution of the finite-time Lyapunov exponents. 

PACS numbers: 05.45.-a,05.45.Pq,05.45.Ac,05.45.Df,95.10.Fh 



I. INTRODUCTION 

Strange nonchaotic attractors (SNAs) are considered 
as typical structures of quasiperiodically forced nonlin- 
ear systems. They are geometrically strange (that is 
they are fractal in nature) just like the chaotic attrac- 
tors, while all their Lyapunov exponents are either zero 
or negative which ensure that the underlying dynamics 
is nonchaotic. Further, due to their fractal nature, the 
SNAs are characterized by aperiodic oscillations. Follow- 
ing the pioneering work of Grebogi et al. [l, SNAs have 
been extensively investigated theoreti cally in several dy- 
namical systems [2^, 4, 5| S 0, 1 i, [Tfl [m, [il [S Q 

m, Ei, [13, m, The existence of SNAs has also 

been demonstrated experimentally [2l|, IH, [l^l in a 
few physically relevant situations. As a consequence, sev- 
eral routes (scenarios having distinct signatures) to SNAs 
have been reported theoretically. These include Heagy- 
Hammel route [l^ , gradual fractalization route [l^ , vari- 
ous types of intermittency routes [l^, [H, [H, H^] , blowout 
bifurcation route @], etc. As all these bifurcation sce- 
narios (routes to SNAs) have been well established in 
the literature, we summarize the different scenarios for 
the formation of SNAs along with their distinct signa- 
tures/mechanisms in Table L Reviews on SNAs can be 
found in Refs. BIH,!!^. 

As mentioned above, while extensive numerical studies 
on the birth of SNAs via different routes are available in 



the literature 

m, [13, m, E 

of them exist 
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ii,iBii[3,[i3,[i3,[il,[i3,[il, 

only a few experimental realizations 
23, 0, [ill- In particular, these exotic 
attractors were confirmed by an experiment consisting of 
a quasiperiodically forced, buckled, magneto-elastic rib- 
bon j23|. SNAs were also realized in analog simulations 
of a multistable potential [l^l, and in a neon glow dis- 
charge experiment [1^ . These attractors were also shown 
to be related to the Anderson localization phenomenon 
in the Schrodinger equation with a quasiperiodic poten- 
tial [13, [1^ . Very recently SNAs have also been observed 
in an excitable chemical system, namely a three elec- 
trode electrochemical cell [3(?|. In this connection, from 
an experimental point of view, nonlinear electronic cir- 
cuits with suitable quasiperiodic forces turn out to be 
especially useful dynamical systems for the identification 
and study of SNAs. For example, Type-I intermittency 
route to SNA was reported in a quasiperiodically forced 
Murali-Lakshmanan-Chua circuit [22!|. Recently, three 
prominent routes, namely Heagy-Hammel, fractalization 
and type-Ill intermittency routes to SNAs, have been 
identified and reported in a quasiperiodically forced nega- 
tive conductance series LCR circuit with a diode [IJ both 
experimentally and numerically by some of the present 
authors. 

In almost all the above studies, as a general rule, the 
driving forces are assumed to be sinusoidal in nature. 
Naturally the question arises as to what happens to the 
dynamics when one or both of the driving forces are non- 
sinusoidal but periodic. Can new routes to the birth 
of SNAs emerge in such a scenario? In order to an- 
swer these questions, we consider the quasiperiodically 
driven negative conductance series LCR circuit with a 
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TABLE I: Routes and mechanisms for the formation of SNAs 

Type of route Mechanism 

Heagy-Hammel [l6{ Collision of period-doubled torus 

with its unstable parent 

Gradual Fractilization [13 ] Increased wrinkling of torus 

without any interaction with 
nearby periodic orbits 

Type-I intermittency [lol ] Due to saddle-node bifurcation, 

a torus is replaced by SNA 

Type-Ill intermittency [S^l Subharmonic instability 

Crisis-induced intermittency Doubling of destroyed torus 
[l^ involves a kind of sudden 

widening of the attractor 
Homoclinic collision [l7| Homoclinic collisions of the 

quasiperiodic orbits 
Blowout bifurcation Due to changes in sign of the 

Lyapunov exponent At transverse 

to the invariant subspace S 
Quasiperiodic route [3, [ij Collision between a stable 

and unstable torus 



diode (which was investigated in [24|) and unravel the 
dynamics of the circuit with one of the forces taken as a 
square wave force (nonsinusoidal) for suitable parameter 
values. The main reason for choosing square wave as one 
of the driving forces is its bistable nature. Bistability is 
responsible for h yste resis in many physical and technical 
systems [sil l32l.l33l Isl j. Further, the square wave has 
also been used for inducing chaos in certain dynamical 
systems [s^. For example a 10 MHz square wave op- 
tical message was injected into a ring laser to produce 
high-dimensional chaotic light 36,]. Thus the study of 
the present circuit has considerable relevance in under- 
standing SNA transitions. 

In the present proposed circuit with a square wave 
force as one of the quasiperiodic forces in addition to 
a sinusoidal force, we have identified a new route for the 
formation of SNA which we term as the bubbling route 
to SNA. In this route bubbles appear in the strands of the 
torus as a function of the control parameter, then the 
sizes of the bubbles increase with the value of the control 
parameter and subsequently strands of the bubbles are in- 
creasingly wrinkled resulting in the birth of SNA ( while 
the remaining parts of the strands of the torus outside the 
bubbles remain largely unaffected) . The mechanism for 
this route is that the quasiperiodic orbit becomes increas- 
ingly unstable in its transverse direction as a function 
of the control parameter which is induced by the square 
wave type quasiperiodic force resulting in an increase in 
the size of the doubled strands (bubbles) in certain parts 
of the main strand and then the doubled strands become 
extremely wrinkled (without a complete doubling of the 
entire main strand) resulting in the SNA. In addition to 
this we have also observed four other prominent routes in 
the same circuit which include fractalization, fractaliza- 
tion followed by intermittency, intermittency and Heagy- 



Hammel routes, the details of which will be published 
elsewhere. 

In order to confirm the existence of the bubbling route 
to SNA in the proposed circuit, we first present a detailed 
numerical analysis of the dynamical equations of the cir- 
cuit in a rescaled form for suitable values of the parame- 
ters using various qualitative and quantitative measures 
to establish this route. These include Poincare surface 
of section, Fourier spectrum, largest Lyapunov exponent 
and its variance, spectral distribution function and dis- 
tribution of finite time Lyapunov exponents. A short ac- 
count of these measures is given in Appendix]^ Next, we 
confirm the results experimentally by the phase portraits 
of the quasiperiodic attractors and SNAs for the corre- 
sponding values of the circuit parameters, again with ap- 
propriate quantification measures, to establish the exis- 
tence of torus and the birth of SNA through the bubbling 
route. 

The paper is organized as follows. In Sec. II, we dis- 
cuss the circuit realization of the quasiperiodically forced 
negative conductance series LCR circuit with diode using 
a sinusoidal and a nonsinusoidal (square wave) forcing as 
quasiperiodic forces. In Sec. HI, we describe the phase 
diagram of the circuit where the regions corresponding 
to the different dynamical transitions to SNAs are de- 
lineated as a function of the control parameters, based 
on our numerical analysis. The birth of SNA via the 
bubbling route as confirmed in the numerical analysis 
is discussed in Sec. IV. Experimental confirmation of the 
bubbling route to SNA is presented in Sec. V. Finally, we 
summarize our results in Sec. VI. Appendix \K\ contains a 
short summary on the identification and characterization 
of SNA and the associated routes. 



II. CIRCUIT REALIZATION 

In this section details about the proposed circuit are 
presented and the circuit equations are written in terms 
of the circuit variables. Then the circuit equations are 
transformed into dimensionless equations (normalized 
equations) using appropriate rescaled variables for a con- 
venient numerical analysis. 

A. Experimental realization: Circuit Equations 

We consider the simple second-order nonlinear dissi- 
pative nonautonomous negative conductance series LCR 
circuit with a sinusoidal voltage generator, fi{t), intro- 
duced by us very recently [24 , |37[ along with a second 
nonsinusoidal force, f2{t), as shown in Fig.[T^. 

The circuit consists of a series LCR network, forced 
by a sinusoidal voltage generator, fi(t), and a nonsinu- 
soidal (square wave) voltage generator, f2{t) (HP 33120A 
series). Two extra components, namely a p-n junction 
diode (D) and a linear negative conductor g^, are con- 
nected in parallel to the forced series LCR circuit. The 
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FIG. 1: (a) Circuit realization of a simple nonautonomous 
circuit. Here D is the p-n junction diode and qn is negative 
conductance. The external emfs are fi{t) = Efi sin(aj/if) and 
/2(i) = Ef2sgn{sm(u}f2t)). The values of the circuit elements 
are fixed as L = 50.3 mH, C = 10.35 nF, R = WOOohms, 
Ef2 = 400my and L0f2 = 17033-ffz. The forcing amplitude 
Efi and its frequency ufi are chosen as the control param- 
eters, (b) i — V characteristics of the p-n junction diode and 
(c) two segment piecewise-linear function. 



negative conductor used here is a standard op-amp based 
negative impedance converter (NIC) . The diode operates 
as a nonhnear conductance, limiting the amplitude of 
the oscillator. In Fig. [1^, v, and in denote the voltage 
across the capacitor C, the current through the inductor 
L and the current through the diode D, respectively. The 
actual v—i characteristic of the diode (Fig.[Tj3) is approxi- 
mated by the usual two segment piecewise-linear function 
( Fig. [T]:) which facilitates numerical analysis consider- 
ably. The state equations governing the presently pro- 
posed circuit (Fig. [J) are a set of two first-order nonau- 
tonomous differential equations 



C— ^iL - Id 
at 



gNV, 



L—— = - KiL - V 
at 



Efi sm{ujfit) 
+ Ef2sgn{sin{ujf2t)), 



where 



ioiv) 



goiv-V), v>V, 
0, v<V. 



(la) 



(lb) 



(Ic) 



Here gjj is the slope of the characteristic curve of 
the diode, Efi and Ef2 are the amplitudes, w/i and 



uj f2 are the angular frequencies of the forcing functions 
fi{t) = Efism{ujfit) and f2{t) = Ef2sgn{sm{ujf2t)), re- 
spectively. In the absence of f2{t), the circuit (Fig. [T^) 
has been shown to exhibit chaos and also strong chaos 
not only through the familiar period-doubling route but 
also via torus breakdown followed by period-doubling bi- 
furcations Here our aim is to investigate the effect 
of the second square wave type external forcing on the 
dynamics and to identify different types of transitions to 
SNAs. 

In order to select the appropriate set of experimen- 
tal parameters for which SNAs can be actually observed, 
we first carry out a detailed numerical simulation (as 
pointed out below) which then serves as a guide and a 
characterizer for experimental investigation. Using such 
an analysis, the values of the diode conductance go, neg- 
ative conductance g^ and break voltage V are fixed as 
1313/iS', — 0.457715* and 0.5^, respectively. We have fixed 
the actual experimental values of the resistance R, in- 
ductance L, capacitance C, external frequency a;/2 and 
forcing strength Ef2 of the square wave as 1900ohms, 
50.3mH, lOMnF, 17033Hz and iOOmV, respectively, 
while we vary the amplitude Efi and the frequency w/i of 
the sinusoidal force as control parameters in order to ob- 
serve the various dynamical states. The forcing functions 
fi{t) and f2{t) are obtained from two separate function 
generators of the type HP33120A. 



B. Numerical analysis: Normalized equations 

In order to study the dynamics of the circuit in de- 
tail, Eq. ([T]) can be converted into a convenient nor- 
malized form for numerical analysis by using the fol- 
lowing rescaled variables and parameters, t = t/^/LC, 
X = v/V, y = (zL/y)(/L7C), El = Efi/V, E2 - 

a = Ry/CjL, 



Ef2/V, 

b = gN\/L/C, and go^/ L/C. 

The normalized evolution equation so obtained from 
Eq. (O is 



= V + f{x). (2a) 

— —x — ay + EiS\i\{4>)+E2sgn{s\i\{9)), (2b) 

= wi, (2c) 

= W2, (2d) 



where 



(h — c)x + c, X > 1, 
hx, X <1. 



(2e) 



Here dot stands for differentiation with respect to t. 
Eq. is then numerically integrated using Runge-Kutta 
fourth order routine to identify the different dynamical 
scenarios corresponding to different values of the rescaled 
parameters. Various interesting dynamical transitions of 
the Eq. (21) are described below. 
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FIG. 2: (Color online) Numerical phase diagram in the 
{El — LUi) plane for the circuit given in Fig.[T] represented by 
Eqs. 3T correspond to period-3 torus, F, BUB, PINT, 
INT, and HH denote the birth of SNAs through gradual 
fractalization, bubbling, fractalization followed by intermit- 
tency, intermittency and Heagy-Hammel routes, respectively. 
C corresponds to the chaotic attractor. Arrows indicate the 
transition regions between two different types of routes to 
SNA. 

III. TWO PARAMETER PHASE DIAGRAM 

The parameter space of the amplitude of the external 
forcing Ei and the frequency loi of the sinusoidal forcing 
is scanned first numerically in the range of Ei E (0.4, 1.1) 
and uji g (1.075, 1.275) to pinpoint different dynami- 
cal behaviors and more specifically for the occurrence of 
SNAs through different routes. From this analysis, var- 
ious dynamical transitions are determined as a function 
of the amplitude of the external forcing Ei and its fre- 
quency LOi. Further, these dynamical behaviors and their 
transitions are also confirmed experimentally for the cor- 
responding values of the experimental parameters of the 
circuit given in Fig. [1] 

A. Numerical analysis 

To start with, we first demarcate the parameter 
space {Ei,uji), by numerically integrating Eqs. into 
quasiperiodic, strange nonchaotic and chaotic regimes 
by using the various qualitative and quantitative mea- 
sures as discussed in the Appendix [X] The numerical 
phase diagram is shown in Fig. [2] for Ei G (0.4, 1.1) and 
u!i e (1.075, 1.275). The various dynamical behaviors in- 
dicated in the phase diagram (Fig. ^ and the interesting 
dynamical transitions are elucidated in the following. 

Transitions from quasiperiodic attractor to SNA and 
subsequently to chaotic attractor occur on increasing the 
value of the amplitude of the sinusoidal force Ei for fixed 
value of its frequency uji . Strange nonchaotic attractors 
created through different mechanisms are found to occur 
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FIG. 3: (Color online) Blow up picture of Fig. [2] in the bub- 
bling transition regime indicated as BUB. 



for different values of the frequency loi of the sinusoidal 
force. Now, we will outline the ranges of values of the 
frequency lui for which SNAs arise from quasiperiodic at- 
tractors through different mechanisms on increasing the 
value of the amplitude Ei of the sinusoidal force. 

Strange nonchaotic attractor created through the 
newly proposed route, namely the bubbling route, is iden- 
tified in the range of frequency uji g (1.085, 1.086). Here 
bubbles appear in the strands of period-3 torus and then 
the bubbles get increasingly wrinkled in the range of the 
amplitude Ei e (0.54, 0.55) of the sinusiodal forcing re- 
sulting in SNA. This phenomenon is named as the bub- 
bling transition to SNA and it is denoted as BUB in 
Fig. [21 A blow up of the two parameter space correspond- 
ing to the bubbling transition is shown in Fig. [31 Further 
increase in the value of Ei ends up in the chaotic behav- 
ior indicated as C in Figs. [21 and [3l Strange nonchaotic 
attractor created through gradual fractalization (F) of 
period-3 (3T) quasiperiodic attractor is identified for 
uji g (1.086, 1.111) as a function of the amplitude Ei E 
(0.5, 0.55). Intermittency route (INT) is found to be ex- 
hibited in the range of frequency uji E (1.111, 1.1268) on 
increasing Ei in the range Ei E (0.5,0.55) and also for 
LUi E (1.1512,1.2615) when Ei E (0.5,0.8). When the 
frequency uji E (1.1268, 1.1512), gradual fractalization is 
followed by intermittency phenomenon on increasing the 
value of the amplitude of the external forcing Ei . It is 
marked as (FINT) in Fig. [21 Torus doubling bifurcation 
from a period-3 torus (3T) to a period-6 (6T) torus and 
then to SNA via the Heagy-Hammel (HH) mechanism 
is found to occur in the range of oji E (1.2501,1.2615) 
on decreasing Ei in the range Ei E (1.1,0.8). The tran- 
sition regions between the above mentioned dynamical 
regimes are indicated by arrows in Fig. [21 which are fixed 
by scanning the frequency uji of the sinusoidal force at its 
fourth decimal place. However, we do not draw a distinct 
boundary between any two scenarios because it requires 
a much detailed numerical analysis on a finer parameter 



scale. 



B. Experimental investigation 

It has additionally been confirmed that the above dy- 
namical behaviors are also exhibited by the experimental 
circuit for the corresponding values of the circuit param- 
eters ^ • and ijJfi{— \J {C/L) ■ uji) by examin- 
ing the two dimensional projections of the corresponding 
attractors obtained by measuring the voltage v across 
the capacitor C and the current il through the induc- 
tor L which are connected to the X and Y channels of 
an oscilloscope, respectively. Here V is the break voltage. 
Then, a live picture of the corresponding power spectrum 
obtained from a digital storage oscilloscope (HP 54600 
series) of the projected attractor has also been used to 
distinguish the different attractors. In addition to this, 
the experimental data of the corresponding attractors 
recorded using a 16-bit data acquisition system [AD12- 
16U(PCI)EH] at the samphng rate of 200 kHz have been 
analyzed quantitatively using the different quantification 
measures, namely the spectral distribution function and 
the distribution of finite time Lyapunov exponents. This 
information is then utilized (i) to pinpoint the different 
dynamical behaviors, (ii) to distinguish the SNAs cre- 
ated through different mechanisms and (iii) also to com- 
pare them with the results of numerical simulation. In 
the following, we will describe only the novel bubbling 
transition in detail by both numerical simulation and ex- 
perimental realization, while the results of other known 
routes will be published elsewhere. 



IV. BUBBLING ROUTE TO SNA: NUMERICAL 
ANALYSIS 

As noted above, in this new route, the bubbles appear 
in the strands of the torus as the value of the amplitude 
El of the sinusoidal forcing is increased for a fixed value 
of its frequency oji . The sizes of the bubbles increase fur- 
ther on increasing the amplitude Ei and the bubbles in- 
creasingly get wrinkled (while the remaining parts of the 
strands of the torus outside the bubbles remain largely 
unaffected) resulting in the birth of SNA. This bubbling 
route is observed in the rather narrow range of frequency 
wi G (1.085, 1.086) as a function of the amplitude of the 
sinusoidal forcing Ei £ (0.5,0.55) indicated as BUB 
in Figs. and [31 It is to be noted that this route is 
significantly different from the well known fractalization 
route [l3| , where the entire strands of the n-period torus 
will continuously deform and get extremely wrinkled as 
a function of the control parameter. The formation of 
SNA through this novel bubbling route has been identi- 
fied in the literature for the first time to the best of our 
knowledge. We have used both qualitative and quan- 
titative measures, which are indicated in the Appendix 
\X[ to confirm the new route. The qualitative proof is 
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FIG. 4: Projection of the numerically simulated Poincare sur- 
face of section of the attractors of Eqs. ((2]) in the (0, x) plane 
for a fixed value of the frequency of the sinusoidal forcing, 
uji — 1.0852, as a function of its amplitude Ei indicating the 
transition from quasiperiodic attractor to SNA through bub- 
bhng route: (a) period-3 torus (3T) for Ei — 0.5, (b) bubbled 
strands of period-3 torus (3T) for Ei = 0.52, (c) enlarged bub- 
bles in the strands of period-3 torus (3T) for Ei = 0.54, (d) 
fractalized bubbles for Ei = 0.546 with the remaining parts 
(away from bubble) of the strands unaffected and (e) chaotic 
attractor (widely interspersed bubbles) for Ei — 0.56. 



given through the Poincare surface of section by distin- 
guishing geometrically between quasiperiodic attractors 
and SNAs. The quantitative confirmation is provided 
using three different measures: (i) The largest Lyapunov 
exponents and its variance are used to distinguish be- 
tween torus and SNA, and SNA and chaotic attractors. 
(ii) Scaling laws deduced from the distribution function 
for quasiperiodic attractors and SNAs are used to dis- 
tinguish them, (iii) Finally, different routes to SNAs are 
also distinguished by the different distributions of local 
Lyapunov exponents. More information on the charac- 
terization is given in the Appendix [X] In the following 
we provide details of the confirmation of the bubbling 
route. 



A. Poincare surface of section plots and power 
spectra 

We have fixed the value of the frequency of the sinu- 
soidal forcing as uji = 1.0852 for illustration and var- 
ied its amplitude in the range Ei G (0.5, 0.55) to elu- 
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(bubbles), as seen in Figs. SJa and [5)3, in certain parts of 
the main strands. This instabihty of the quasiperiodic at- 
tractor arises due to the presence of the square wave pulse 
(finite amplitude for finite durations). Further increase 
in the value of the amplitude of the forcing {E\) results 
in an increase in the size of the doubled strands (bub- 
bles) as shown in Figs.lH^ and[SJ;, and then the doubled 
strands become extremely wrinkled (without a complete 
doubling of the entire main strand) resulting in the SNA 
as depicted in Figs. [Hi and[SJl. 

We now provide quantitative confirmation of the above 
results to distinguish between torus and SNA, and SNA 
and chaos. 
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FIG. 5: Enlarged figures of Figs. |4]to show the bubbling tran- 
sition to strange nonchaotic attractor. 



cidate the emergence of bubbling route to SNA in the 
present system ([2]). The Poincare surface of section plot 
of the three strands corresponding to period-3 torus for 
the value oi Ei = 0.5 is shown in Figs. [4^ and[5K- The 
corresponding phase portrait and power spectrum are de- 
picted in Figs.[6^(i) and[6^(ii), respectively. As the value 
of the amplitude E\ is increased further, bubbles start to 
appear in all the three strands starting from E\ = 0.516. 
These are shown in Figs. [3Jd and [3)3 for E\ = 0.52 and 
the corresponding phase portrait and power spectrum are 
shown in Figs. [5)3(1) and[5)3(ii), respectively. Further in- 
crease in the value of Ei results in an increase in the size 
of the bubbles as shown in Figs. [4); and [5); for the value 
of El = 0.54, whose phase portrait and power spectrum 
are shown in Figs.[6l:(i) and[6j;(ii), respectively. Beyond 
the value of Ei = 0.54, the strands of bubbles deform 
and get increasingly wrinkled (while the other parts of 
the strands of period-3 torus outside the bubbles remain 
unaltered as seen in Fig. [SJl) leading to the formation of 
SNA as depicted in Fig. [3)1 for the value of Ei 0.546. 
The phase portrait and power spectrum for this value 
of El are shown in Figs. [B)i(i) and[5)l(ii), respectively. 
Finally, to confirm that the SNA transits to a chaotic at- 
tractor beyond Ei = 0.55, we have depicted the Poincare 
surface of section of the latter in Figs. [4^ and[5]3 with the 
corresponding attractor and power spectrum in Figs. [6^ 
for El = 0.56. 

The mechanism for the bubbling route is that the 
quasiperiodic orbit becomes increasingly unstable in its 
transverse direction as a function of the control parame- 
ter (El ) , resulting in the formation of the doubled strands 



B. Largest Lyapunov exponent and its variance 

The largest Lyapunov exponent. A, and its variance, 
/i, that is the variance of A from finite time Lyapunov 
exponents Ai(7V)'s, i — 1, 2, • • • , M of length N, defined 
as 



M 



(3) 



are shown in Figs. [7] in the range of Ei £ (0.54,0.546). 
The attractor depicted in Fig. [6)l(i) for the value of 
El — 0.546 is strange but it is nonchaotic as evidenced 
by the negative value of the Lyapunov exponent shown in 
Fig. [7^. It is also to be noted that both the Lyapunov ex- 
ponents and its variance (Fig. [T)^) clearly indicate a crit- 
ical value of ampfitude £'f = 0.5432, {Ei < Ei), below 
which torus exists and above which, {Ei > Ei), SNA ap- 
pears. The regions of torus and SNA are clearly indicated 
by smooth and irregular variations, respectively, in the 
values of both the Lyapunov exponents A and its variance 
fi. Finally the transition of SNA into a chaotic attractor 
is confirmed by the change in the largest Lyapunov ex- 
ponent from negative to positive values at Ei = 0.55 as 
shown in the inset of Fig. [7^. 



C. Spectral distribution function and scaling laws 

In order to distinguish further whether the attractors 
depicted in Figs. [5] are quasiperiodic or strange non- 
chaotic or chaotic attractors, we proceed to quantify the 
changes in their power spectra. The spectral distribution 
function, defined as the number of peaks in the Fourier 
amplitude spectrum larger than some value cr, is used to 
distinguish between quasiperiodic attractors and SNA as 
well as SNAs and chaotic attractors. The quasiperiodic 
attractors obey a scaling relationship N{a) ~ \ogiQ{l/a), 
while the SNAs satisfy a scaling power-law relationship 
N{a) - cr~'^,l < P < 2 ^. Similarly for the chaotic 
attractor, the scaling relation is N{<t) ^ <t^'^,(3 > 2. 
Spectral distribution functions (filled circles) of the torus 
(Fig. [5^) and bubbled torus (Fig. [B))) satisfy the scaling 
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FIG. 6: Projection of the numerically simulated attractors 
and their power spectrum of Eqs. Q for the same values 
of the frequency u-i and the amplitude E\ of the sinusoidal 
forcing as in Figs. ID (a) period-3 torus (3T), (b) bubbled 
period-3 torus, (c) period-3 torus with enlarged bubbles, (d) 
fractalized bubbles (SNA) and (e) chaotic attractor : (i) phase 
portrait in the (x, y) space; (ii) power spectrum. 



relation N{a) ^ log]^Q(l/cr) as indicated by the solid line 
in Figs. [5^ and[S)3, respectively, which is the characteris- 
tic of a torus. On the other hand the spectral distribution 
function of the SNA (Fig. exhibits power-law behav- 
ior as depicted in Fig. [5J; (filled circles) with the value of 
the exponent /3 = 1.88, confirming the existence of SNA. 
For the chaotic attractor (Fig. [6^), the scaling exponent 
(Fig. [SJi) turns out to be /3 = 3.5 as required. Again the 
solid lines in Figs.[Ht and[5Jl represent the scaling law for 
SNA and chaos, respectively. 




0.54 0.542 0.544 0.546 

El 




0.54 



0.542 0.544 0.546 



FIG. 7: Transition from torus to SNA through bubbling 
route for the same value of frequency as in Figs. |4] and in 
the range of amplitude E\ £ (0.54, 0.546) obtained numeri- 
cally, (a) Largest Lyapunov exponent (A) and (b) its vari- 
ance (fj,). Inset in (a) depicts transition from SNA to chaos 
for El G (0.545,0.56). 
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FIG. 8: Spectral distribution function (filled circles) calcu- 
lated numerically, (a) torus (Fig. [S^), (b) bubbled torus 
(Fig. ISJd), (c) SNA (Fig. |6ll) and (d) chaotic attractor 
(Fig. |6t). Solid lines in (a) and (b) corresponds to the scal- 
ing relation N{a) ~ logj^Q(l/cr) and in (c) and (d) correspond 
to the scaling relation N(a) ~ o-^*^, with /3 = 1.88 and 3.5, 
respectively. 
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FIG. 9: Distribution of finite time Lyapunov exponent calcu- 
lated from both numerical data (solid line) and experimental 
data (dashed line) of (a) torus (Fig.Igli) and (b) SNA (Fig.[li). 



V. BUBBLING ROUTE TO SNA: 
EXPERIMENTAL CONFIRMATION 

As a next step, in order to confirm the results of our 
numerical simulation in the experimental circuit shown 
in Fig. [1] a snapshot of the dynamical behavior for the 
corresponding values of the experimental parameters is 
obtained (as mentioned in Sec. |TT| and compared with 
that of the numerical results. Further, the correspond- 
ing experimental data are analyzed using various quan- 
tification measures mentioned in the previous section to 
confirm the nature of the dynamical behavior. 



D. Distribution of local Lyapunov exponents 



In addition to the qualitative discussion through the 
Poincare surface of section plots in the (0, x) plane 
(Figs.mandO in distinguishing the type of route through 
which SNA appears, it is also possible to distinguish 
the same using the distribution of a quantitative mea- 
sure, namely finite time Lyapunov exponents. It has 
been shown that a typical trajectory on a SNA 

actually possesses positive Lyapunov exponents in finite 
time intervals, although the asymptotic exponent is neg- 
ative. As a consequence, it is possible to observe different 
characteristics of SNAs created through different mech- 
anisms by a study of the differences in the distribution 
of finite time exponents P{N,X) The distribution 

can be obtained by taking a long trajectory and divid- 
ing it into segments of length TV, from which the local 
Lyapunov exponent can be calculated. In the limit of 
large N, this distribution will collapse to a 5 function 
P{N, A) — > 6{A — A). The deviations from— and the ap- 
proach to— the limit can be very different for SNAs cre- 
ated through different mechanisms p^ . 

We have calculated the distribution of local Lyapunov 
exponents P{N,X), for N = 2000, for the attractors 
shown in Figs. [n^(i) and EJi(i) in order to confirm the 
nature of transition to SNA. The distribution of the local 
Lyapunov exponents for the period-3 torus (solid line) is 
shown in Fig. [9^ in which the local Lyapunov exponents 
are peaked about the largest Lyapunov exponents (neg- 
ative values) of the torus while that of the SNA shown 
in Fig. [Hb by solid line has its maximum at a positive 
value of the local Lyapunov exponents. The distribution 
of local Lyapunov exponents for SNA exhibits an elon- 
gated tail in its negative values because of the fact that 
in the bubbling transition parts of the strands of period- 
3 torus remain unaffected even after the birth of SNA 
which contributes largely to the negative values. This 
confirms the existence of bubbling transition to strange 
nonchaotic attractor. 



A. Phase portraits and power spectra 

We have depicted the snapshots of the phase portraits 
and the corresponding power spectra of the attractors as 
seen in the oscilloscope (which is connected to the cir- 
cuit shown in Fig. [T]) in Fig. [10] for the corresponding 
values of the parameters of numerical simulation. Ex- 
perimental period-3 torus and its power spectrum corre- 
sponding to the numerical results. Figs. [5^, are shown in 
Figs. fTOb(i) and [TOb(ii). The attractors in the bubbling 
regime for the values of the amplitude of the sinusoidal 
forcing Eji = 0.26V and 0.27V are shown in Figs.[TUb(i) 
and c(i), respectively. The corresponding power spectra 
are shown in Figs. fTOb (ii) and c(ii), respectively. Experi- 
mental phase portrait of the strange nonchaotic attractor 
and its power spectrum for the value of Efi = 0.273^ 
are depicted in Figs.fTUHfi) andfTUtlfii). respectively. It is 
also seen that the spectra of the quasiperiodic attractors 
are concentrated at a small discrete set of frequencies 
while the spectrum of the SNA has a much richer set 
of harmonics. Further the resemblance of the attractors 
illustrated in Figs. \TO\ i] with that of the attractors in 
Figs-in^i) confirms the existence of bubbling transition to 
SNA in this negative conductance series LCR circuit with 
diode having both the sinusoidal and nonsinusoidal forces 
as quasiperiodic forcings. Finally, the chaotic attractor 
for Efi — 0.28V and its power spectrum are shown in 
Figs. [Tl3. 



B. Spectral distribution function and scaling laws 

In order to confirm that the experimental phase por- 
traits shown in Figs. [TUb . [TUb . [TUt l and \Wk are indeed 
that of torus, bubbled torus, SNA and chaotic attractor, 
respectively, the corresponding data are examined for the 
behavior in their spectral distribution. Figs. [Tib andlllb 
show the spectral distribution function (filled triangles) 
for the torus in Figs. [TUb and [TUb satisfying the scahng 
relation N{a) ^ logj^g (!/''') indicated by the solid lines 
while that of the SNA (Fig. [TOli) shown in Fig. [TTt obey 
power-law distribution with the value of the exponent 
(3 = 1.96 lying within the characteristic value for SNAs. 
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FIG. 10: Snapshots of the experimental attractors and their 
power spectrum of the circuit shown in Fig. [1] for the corre- 
sponding values of the frequency u/i and the amplitude 
of the sinusoidal forcing in Fig.|4l (a) period-3 torus (3T), (b) 
bubbled period-3 torus, (c) period-3 torus with enlarged bub- 
bles, (d) fractalized bubbles (SNA) and (e) chaotic attractor: 
(i) phase portrait in the {vcil) space; (ii) power spectrum. 



For the chaotic attractor (Fig.fTOb). the scahng exponent 
turns out to he P — 4.0 (Fig. fTTH) as expected. 



C. Local Lyapunov exponents 

Further, in order to examine whether the SNA shown 
in Fig.fTOH arises from the bubbhng transition, the distri- 
bution of the local Lyapunov exponents calculated from 
the experimental data of the torus (Fig. flOk ) and the 
SNA (Fig. [inii) are depicted in Figs. [SJi and [Hh respec- 
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FIG. 11: Spectral distribution function (filled triangles) cal- 
culated from the experimental data of (a) torus (Fig. llOh 'l. 
(b) bubbled torus (Fig. [lOb), (c) SNA (Fig. [lOji) and (d) 
chaotic attractor (Fig. llOb ). Solid curve/line in (a) and (b) 
correspond to the scaling relationship for the quasiperiodic at- 
tractors and in (c) and (d) correspond to the scaling relation 
for the SNA and chaotic attractor, respectively. 



tively as dashed hues. The elongated tail in the distri- 
bution of the local Lyapunov exponents even for SNA 
(Fig. [TUH-) confirms the existence of undisturbed strands 
as shown in Fig. thereby confirming experimentally 
the birth of SNA via the bubbUng transition. 



VI. SUMMARY AND CONCLUSION 

In this paper, we have reported the birth of strange 
nonchaotic attractors through a novel route which we 
term as the bubbling route to SNA in a negative con- 
ductance series LCR circuit with the diode containing 
nonsinusoidal (square wave) force as one of the quasiperi- 
odic forcings. At first, we have presented the numerical 
analysis of the dynamical system, namely Eq. ^ of the 
circuit (Fig. [1]) for suitable ranges of the amplitude Ei 
and the frequency ui of the sinusoidal force while the 
other parameters are held fixed. Following this, we have 
also confirmed the numerical results experimentally by 
the snapshots of the phase portraits of the quasiperi- 
odic attractors and SNAs as well as chaotic attractors 
for the corresponding values of the circuit parameters. 
Further, the numerical and experimental data have been 
analyzed using various quantification measures attribut- 
ing to the existence of torus, SNA, birth of SNA through 
the bubbling route and transition to chaos. In particu- 
lar, we have characterized the quasiperiodic attractors, 
SNAs and chaotic attractors using maximal Lyapunov 
exponent and its variance, Poincare maps, Fourier am- 
plitude spectra, spectral distribution function and distri- 
bution of finite time Lyapunov exponents. The distri- 
bution of local Lyapunov exponents indeed clearly dis- 
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tinguishes the characteristic properties of both the torus 
and the SNA, confirniing the existence of bubbUng route 
to the SNA. The experimental observations, numerical 
simulations and characteristic analysis showed that the 
simple dissipative negative conductance series LCR cir- 
cuit even with a nonsinusoidal (square wave) force as 
one of the quasiperiodic forces does indeed admit strange 
nonchaotic behaviors of different types and in particular 
admits a novel bubbling route to SNA. 
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APPENDIX A: IDENTIFICATION AND 
CHARACTERIZATION OF SNAS AND THEIR 
ROUTES 

Torus, SNA and chaotic attractors and the transi- 
tions between them through different routes can be 
identified and characterized through various qualitative 
and quantitative measures. In this Appendix, we sum- 
marize the main measures used in the recent litera- 
ture 

@, i, [M Q E, d, US, [U m, [13, m in the analy- 
sis of transitions to SNAs from torus attractors and from 
SNAs to chaotic attractors. In the present work also, we 
utilize these measures. 

1. Qualitative measures: 

Geometrically smooth (torus) and non-smooth 
(SNAs and chaotic) attractors can be distinguished 
qualitatively using Poincare surface of sections and 
Fourier spectra. The Poincare surface of section 
shows smooth strands for quasiperiodic attractors, 
non-smooth strands for SNAs, widely interspersed 
points throughout the phase space for chaotic at- 
tractors, which clearly reveals whether an attractor 
possesses a geometrically smooth or complicated 
structure. The spectra of the quasiperiodic attrac- 
tors are concentrated at a small discrete set of fre- 
quencies while the spectra of SNAs and chaotic at- 
tractors have a much richer set of harmonics. 



Further, different types of routes to SNAs and their 
mechanisms for their formation can also be identi- 
fied qualitatively using the Poincare surface of sec- 
tions by observing the nature of the dynamics in 
these plots as a function of the control parameter. 
Different routes for the formation of SNAs have dif- 
ferent characteristic dynamics in their Poincare sur- 
face of section. 



2. Quantative measures: 

(a) The largest Lyapunov exponents can be used 
to distinguish between (i) torus and SNAs and 
(ii) SNAs and chaotic attractors. Torus mo- 
tion is characterized by a smooth negative 
Lyapunov exponent, SNAs are characterized 
by either zero or non-smooth negative Lya- 
punov exponents as a function of control pa- 
rameters and chaotic attractors have atleast 
one positive Lyapunov exponent. Further, the 
transition from torus to SNAs exhibits differ- 
ent signatures in the values of the largest Lya- 
punov exponents and their variance for differ- 
ent routes to SNAs ^ . 

(b) Further, torus and SNA can also be dis- 
tinguished quantitatively by using the spec- 
tral distribution function, which is defined as 
the number of peaks in the Fourier ampli- 
tude spectrum larger than some value a [3|. 
The quasiperiodic attractors obey a scaling 
relationship N{a) ~ \ogiQ{l/a), while the 
SNAs satisfy a scaling power-law relationship 
N{a) ^ a~f^, 1 < /3 < 2. For chaos, the scal- 
ing exponent /3 > 2. 

(c) Finer distinction between the different types 
of routes for the formation of SNAs can also 
made using the distribution of finite time Lya- 
punov exponents. Different routes are charac- 
terized by different types of the distribution of 
finite time Lyapunov exponents [10]. 



The different signatures of the above quantitative mea- 
sures corresponding to different scenarios (routes) for the 
formation of three well known types of SNAs are tabu- 
lated in Table II. 
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